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ABSTRACT 

A model for a supersonic blade row with two In-passage shock waves Is 
developed. It accounts for three-dimensional effects In observed flows by 
using a hypothetical blade shape In a two-dimensional cascade. There Is 
enough flexibility In the choice of blade shape to accommodate a desired 
entrance angle, exit angle, boundary-layer thickness and stage pressure ratio 
at a given entrance Mach number. 

The model divides the mean flow Into regions of uniform or one-dimensional 
flow In which the solutions for the superimposed unsteady flow may be found 
successively. The analysis makes use of previous solutions for unsteady flow 
In cascades and over an oscillating wedge. 

Six flow conditions are chosen In the range of parameters for which the 
two-shock model Is valid for studies of flutter In torsion and bending. It Is 
found, In keeping with previous results from a single-shock model, that In 
each case there Is Increasing Instability with decreasing frequency. 


INTRODUCTION 

The analysis of unsteady flows In blade rows has been facilitated by the 
use of models for the underlying steady flow which permit the governing differ- 
ential equations to be written with constant coefficients. This has been true 
for both supersonic flows and mixed flows with a single shock wave (ref. 1) In 
which the mean flows were assumed to be uniform. A number of experimental 
measurements (refs. 2 to 5), using both pressure sensors and laser anemometers, 
has shown that multiple shock waves are a common occurrence In the tip region 
of the blades, especially near the operating point of the compressor. One of 
the frequent patterns (fig. 1) observed In the experimental studies has an 
oblique shock wave at the leading edge of each blade and a nearly normal shock 
wave at the trailing edge. It Is that flow which we shall Investigate for 
bending and torsional flutter. 

We shall first Introduce a model for the mean flow which Incorporates the 
two shock waves and modifies the actual blade surfaces to account for three- 
dimensional flow effects. There are three mean-flow domains which are modeled 
as either uniform or one-dimensional flow so that the greatest use can be made 
of previous solutions for unsteady duct flow (ref. 1) and flow around an oscil- 
lating wedge (ref. 6). The range of validity of the model Is estimated and six 
numerical examples are calculated covering the Intervals of Inlet Mach number 
and stage pressure ratio for which the model Is valid. The stability charac- 
teristics of the examples are discussed. 



Considerable use is made of the analysis previously developed in 
reference 1 for the flow in a blade passage with a single normal shock wave 
present. Some familiarity with that paper is assumed in parts of the exposi- 
tion. The stability analysis of a mean-flow configuration related to those 
treated here and in reference 1 has recently been completed (unpublished data 
by E. Acton, Unsteady Flow in a Supersonic Cascade with Strong Oblique Leading 
Edge Shocks (1982)); it treats the case of a single strong shock wave 
emanating from the blade leading edge. It uses the principal elements of the 
blade model employed here. 

A computer program, TIPS (Two In-Passage Shocks), has been written which 
uses the analysis described in this report to calculate the stability of a flow 
In a cascade. The user provides inlet flow conditions, geometric parameters 
and disturbance characteristics and the program yields corresponding unsteady 
lifts and moments. Contact COSMIC, The University of Georgia, Athens, Georgia 
30602, concerning the availablity of this program. Frank Molls programmed the 
model and supplied the numerical calculations. 


MODEL 

An example of a double-shock configuration observed in an experimental 
rotor (ref. 5) is shown in figure 2 which exhibits static pressure contours 
measured at the blade tips. The dashed line beginning at the leading edge of 
a blade corresponds to an oblique shock wave whose angle is calculated from 
the pressure ratio across the contours. The position of the shock is in good 
agreement with the direction of the contours. A nearly normal shock wave ema- 
nates from the trailing edge of the blade. The apparent discrepancy between 
the angle of the oblique shock wave and the small wedge angle of blade at its 
leading edge is attributed in reference 5 to a radial contraction of the flow 
in the blade passage. 

Following the course set in an earlier analysis (ref. 1) we wish to use a 
model of the flow in the blade passage which represents the main features of 
the mean flow to a good approximation and at the same time keeps the equations 
governing the unsteady motion as amenable to solution as possible. In figure 
3 the incremental rotor annulus at the blade tips has been unrolled into a 
two-dimensional cascade. In it a two-dimensional contraction is Imposed by 
the presence of a wedge at the leading edge. The wedge angle, 6 , while much 
larger than the metal angle of the blade, correctly accounts for the shock 
angle, * . Where the oblique shock intersects the lower blade, the surface of 
the blade is turned parallel to the flow behind the shock. Consequently the 
shock is not reflected, in accordance with the observed flow in figure 2 . 

In order to terminate the blade In a small thickness comparable to the 
combined thicknesses of the upper- and lower-surface boundary layers, the upper 
surface of the blade after its intersection with the oblique shock wave is 
represented by a (cubic) function f u (X 3 ) whose initial direction is that 
of the turned flow, 83 , and whose final slope corresponds to the exit 
angle, 83 , which is a free parameter of the blade design. Likewise, the 
lower surface over the same range of the coordinate X 3 is designated by 
the function f a (x 3 ). The shock wave terminating this region of flow is 
normal to the exit direction, U 2 - The last segment of the lower surface 
along which the flow is subsonic following the second shock wave, is parallel 
to the exit flow direction. 
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The mean flow used In this model Is uniform In the entrance region, the 
triangular region Immediately behind the oblique shock wave and In the exit 
region behind the normal shock wave. Solutions to the unsteady equations of 
motion can, therefore, be obtained analytically In those regions. The flow In 
the channel between the curvilinear surfaces Is assumed to be one-dimensional 
so that the unsteady motion In that region may be obtained subject to a slowly- 
varying flow approximation. 

The unit of length In the model Is the length of the blade as measured In 
the turned flow direction 63. The chord length Is, then, 

l - 1 > 1 

c cos(x 0 - 63) - 


Other lengths are 


1 + ' 
s = 8 cos B./sol, s =8. sin 13 , /sol 
cl cl 

+ 1 

s g = s s 1 n(x - S)/s 1 n x, s 3 = 8 c sin ( 3 3 /sol 

d $ = 1 - (s + + s cot x) cos 6, 8^ = (1 - d $ ) tan 6 

(Here "sol" denotes the solidity of the blade row.) 

The tralllng-edge values of the boundary-layer displacements are specified 
In fractions 6 U and 5 ^ of the base width, 8^, °f a blade with unidirec- 
tional flow after the oblique shock. The chord Intersects the base at a 
fraction v of Its width measured from the lower side. 

sol sin x sin (x Q - ( 3 3 ) 
cos(( 3 1 - x) sin 6 

The two cubic functions defining the curved portions of the blade surface 
are 



The cross-sectional area of the blade passage per unit blade height of the 
cascade Is 
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In these expressions , 



and 


d a = d s ‘ s 2 s1n ( B 2 " |3 3 ) 

A convenient but not necessary assumption In this model Is that the 
elastic axis for torsion lies on the chord at some fraction f of Its length 
measured from the leading edge. Then the distance d 0 of the elastic axis 
from the origin of x-| In the Inlet direction Is 

d = -s + + f[cos 6 + v(l - d ) sin 6 tan 6] 
o s 

which reduces to the value appropriate for a system of Infinitesimally thin 
blades (ref. 1) when 6=0. 

In region 2 the corresponding moment arm of the elastic axis Is 
d o,2 = 1 ~ s 2 “ V " cos ( |3 2 " x o ^ 

When allowance for the tralllng-edge boundary-layer thicknesses Is made, the 
blade separation and stagger lengths are 

I I I 

s 2 = cos 0 2 /sol - (6 u + 6^) a b cos(B p - & 3 ), 

s p = i c sin (3 2 /sol - (6 u + 6 a )i b sin(0 2 - |3 3 ) 

At any entrance Mach number and for any cascade geometry this model 
Is valid If the static pressure ratio P 21 across the blade row Is not so low 
that the oblique shock falls to Intersect the adjacent blade or so high that 
the velocity behind this shock Is subsonic. Figure 4 shows the region of 
validity In the M-| , P 21 plane for a cascade of solidity 1.4, stagger 
angle 62° and entrance-flow angle 67°. 1 To the left of the valid area Is the 
region of low pressure ratio In which the flow through the blade passage Is 


^The blades are untapered at the trailing edge and take the shape Indi- 
cated by the dashed lines In figure 3. 
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likely to be entirely supersonic. To the right Is the region of high pressure 
ratio In which a strong oblique shock Is a likely configuration. 

In general, even though the blades are tapered near the trailing edge to 
more nearly approximate the shape of real blades, there Is a finite tralllng- 
edge thickness, corresponding to a boundary-layer displacement thickness, and 
a wake of finite thickness as well. It Is assumed. In keeping with the bound- 
ary conditions on the Infinitesimally thin wake In reference 1, that the 
pressure and normal component of velocity are also continuous across the wake 
of the present model . 


ANALYSIS 

The unsteady flow produced by the oscillating blade row In the supersonic 
region upstream of the leading shock has been treated In reference 1 for the 
case where there Is a single normal shock at the leading edge of each blade. 

In the present double-shock configuration the upstream supersonic region 
extends an additional length s cot * along the upper blade surface before 
the oblique shock from the adjacent blade touches down. The upstream solution 
Is the same as In reference 1 and extends without modification through this 
region up to the leading shock. Likewise, the unsteady solution In the sub- 
sonic flow region 2 downstream of the normal shock wave Is, with slight modi- 
fication, the same as the downstream subsonic solution of the single-shock 
analysis of reference 1. The main new features of the present model are the 
oblique shock and the supersonic region between this shock and the trailing 
normal shock. 


Boundary Conditions on the Blade 

With the bending motion of amplitude H 0 assumed normal to the entrance- 
region surface and with x = d 0 the projection of the elastic axis upon that 
surface, the amplitude of the upper surface oscillation of the zeroth blade Is 


The motion of the lower surface of the same blade Is 

W 0 = H cos 6 + A (x 0 + 1 - d - f) 
3 o o 3 s 


(Subscript 1 refers to the entrance flow, subscript 3 to the region behind the 
oblique shock. Subscript 2 Is reserved for the subsonic exit region. See 
fig. 3.) 

As Indicated In reference 7, the most general blade motion Is the super- 
position of motions with constant Interblade phase angle, a, so we need only 
consider motions of this type. Then In the cascade, the motion of the n^ 
blade Is determined by 

V 3 ( x 3 * ns 3-»3 * "<‘b * s 3>) - elT "’ V 3<*3' y 3 ) (1) 
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where V 3 can be any of the physical variables. In particular, on the lower 
surface of the first blade, 

W 3 = e io [H 0 cos 6 + A q (x 3 + as 3 - f)] 

(We shall see later that there Is an additional motion with a. = cot(x - <$) 
associated with the perturbation of shock position. However, It does not enter 
Into the shock relations which are the Immediate concern of this discussion.) 


Conditions on the Oblique Shock Wave 

As discussed In reference 1, the pressure disturbances generated In the 
upstream supersonic region distort the oblique shock and may be transmitted by 
It, but not reflected, l.e., the flow region 1 Is unaffected by the shock. 

The perturbed Ranklne-Hugonlot equations have been given by Carrier (ref. 6 ) 
for an oblique shock In a disturbance-free upstream flow, and by Moore (ref. 8 ) 
when there are upstream disturbances. A particular case of the latter Is given 
In reference 1 for the normal shock wave. It also applies to oblique shock 
waves when the coordinates move with the component of the mean velocity paral- 
lel to the shock. The dimensionless time derivative In the moving coordinates 
transforms as 


- cos x 

since time was normalized by c/^i sin x- With the displacement of the shock 
normal from Its mean position denoted by x s (y,t), the pressure perturbation by 
p 1 , the disturbance normal and tangent velocities by u( n ) and u(t), respec- 
tively, then the shock conditions, to first order In the disturbance quan- 
tities, are 


. 1 

ay j ~ sin x Dt 


- — » — - — (— 
at sin x l at 


p 3 = 


(2H* s 1 n 2 x - v + A , 4 [ (n) 

j p i + iT+T s1n x |^ n) - 


Vs 


ot 




Vs 

(v - 1)M 2 s1n 2 x - 2 

L(n) °s x s 


Dt = 

( y + 1 )M 2 sin 2 x 

[ U 1 ■ Dt 

t 2 v. 


-~-l\ _D_ 
\v + 1 / sin 


4 " - u i (t) 


# 3 s 1 n( x - 5) 
^ sin x 


ax 

c 

ay] 


sin x 


where, as Indicated previously, the subscripts 1 and 3 refer to conditions 
upstream and downstream of the shock. When these relations are expressed In 
terms of the streamwlse and cross-stream velocities In the upstream and post 
shock flows 
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u.| = sin x - uj"^ 


cos x, v 


1 


sin x + uj n ^ 


cos x 


u 3 = u< n) s1n(x - 5) - u< l) cos ( x - 6), v 3 = u^ t) s1n(x - 6) + u^ n) 
they take the form 

D s x s 

p 3 + ! 1 ^ "Dt 1 = “ll p l + a 12 u l + a 13 v l 
% D s x s 3x s 

U 3 + l 2d^ Dt + K 2 3y = a 21 P l + a 22 U l + a 23 V l 

o S 

^1 D s x s 3x s 

V 3 + *3 ¥ z ~W + K 3 iy^ = a 31 P l + a 32 U l + tt 33 V l 


cos(x - 6) 

( 2 ) 

(3) 

(4) 


where we have equated tangential velocities at the shock wave to obtain 
cos(x - 6)/cos x =^V%- The coefficients I m , K m , a pq are the functions 
of the entrance Mach number and shock angle listed In appendix A. The quan- 
tities on the right side are known from the solution for the unsteady flow In 
the upstream supersonic region which, as we Indicated, Is Independent of the 
downstream flow. 


Solution In Region 3 

Carrier (ref. 6) developed the solution to the small amplitude, unsteady 
flow behind a shock wave attached to a wedge that Is oscillating about Its apex 
as Bessel series for the unsteady potential, vortlclty and shock displacement. 
We adapt his solution to the region behind the oblique shock by generalizing 
the wedge motion and replacing a portion of the potential series representing 
the effect of the boundary condition on the blade by a form characteristic of 
duct solutions as, for example, In reference 1. The potential, stream func- 
tion, and shock displacement In region 3 between the zeroth blade and the first 
blade (figure 3) are, In our notation, 


1k„M„(x„+as„) 

$ 3 = e 2-i cosh ne J n ( k 3 r ) + 

n=l 


(5) 


'V 


3 


1w 3^ X 3 +aS 3 M 3 /B 3~ ay 3 /l3 3^ 






( 6 ) 


7 



(7) 


1 " 3 “ M 3< s 3- y 3> /B l V'jj 

n n 
n =0 


where 

r = (x 3 + <ss 3 ) - fi 3 ( s 3 - y 3 ) 



e = tanh 


B 3^ S 3 y 3^ 

(x 3 ♦ as 3 ) 


The potential §5 Is given by 

*b ■ s K c ° s s * \< x 3 * 1 - “s - f )]f,(y 3 ;«3» 3 ) 



3co 3 w 3 f V y 3 ;l W + w 3 f 2^ y 3 ;M 3 w 3^ 


co 

CO 


1 

3 


In which the functions f 1 ( y ;co) and f 2 ( Y *“) are those already defined for 
region 1 in reference 1 and represent duct effects having contributions from 
all the blades. We have followed Van Dyke (ref. 9) In adding to the series 
for the shock displacement the term 

VokV * 2 - 0 3 < s 3 - *3>] 

which describes the change In shock position due to the motion of the blade 
leading edge. When the blade has both a plunging motion and an oscillation 
around an axis not at the apex, 

D q = e 1o |n 0 cos x - A 0 fc c f cos(x + x 0 - B 1 )J (8) 

The perturbation pressure and velocities are obtained from the velocity 
potential and stream function from the relations 



From the first of these and the periodicity relation (1) the pressure on the 
blade surface behind the oblique shock wave Is 
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f V X 3' -8 V " 


k 3 s 1[k 3 H 3 (x 3 +l-d s )-o] 
2s 7 e 


Ei 

n=0 


A + — A - A 
H n+1 M„ n rt n-l 

v 3 ) 




X 3„|k 3 (x 3 + 1 - d s ) 


+ A (x~ + 1 - d r + s-j 
O' 3 S3 


)] 

3 * f >] 


-la 
<o, e s 


„ ,co„ H COS 6 

S„ I 3 0 


+ 1A 0 ( f i ( s 3^ M 3 w 3) 


1co,A S . ( M„(i)_s„ | 

s 3 e r3 f 2 (s 3 ;M 3 u 3 ) ~ s1n(M 3 co 3 S 3 ) f l (0;M 3 u 3 ) j 


d r - 1 < x„ < B,s, - s Q 
S 3 3 3 3 


The coefficients A , C , D remaining unspecified In the solutions (5), 

n n n 

(6), (7) must be determined from the shock conditions (2), (3), (4) In a manner 
to be described below. 

For convenience In the later computation of the lift and moment on the 
blade we record here the contributions to those quantities obtained from the 
foregoing expression for the pressure. 

The lift developed on the blade surface from the leading edge to 

x 3 ’ - s 3 1s 


+ 



- aSu. 


aS, 


u 3 \ H 3 + 2 A o) M S 3 ;M 3°V + 1A o H 3^ S 3^ 


-la 


where 
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=. ■ r 


1M 3 k 3 n 


<n+lr^y(2il + n + l,-1M,k,as) 

J_( k 3 n) d n = in- > — 

n 3 ^3 k 3 (2H 3 ) 2J ' fi.1 (n + 8.) ! 


H_ = H cos 5 + A ( s„ - f) 
3 0 0 3 


and 


H 3 (s 3 ) _ f 2 (s 3 :M 3“3 ,S 3 )a3 3 + f 1 (s 3 ;M 3 u 3 ) - 

The evaluation of a n Is made by expanding the Bessel function In Its power 
series at n = 0 and Integrating to obtain a sum of Incomplete gamma functions 
t( 29. + n + 1 ,-1M 3 k 3 as 3 ) . 

Likewise, the moment contributed by this blade surface, measured about the 
leading edge, Is 


p 

M 33 ~ I 

V 1 


(x 3 d $ + 1 )P 3 (x 3 ,-2. b )dx 3 


E (Vl * S; % - A n-l)§ n 


k a s , 

517 Z, ( A nti * 5" 

J n=0 


2 -1«J r- 


ss 3 e 



2 N 

A 

1 A 

[(" 

3 + A o 3 aS 3, 

| o> 3 f -j (s 3 ;M 3 <o 3 ) + 1A 0 H 3 (s 3 ) 


where 


aS, 

■ 


1M„k„n 

ne J J J n (k 3 n)dn 


00 

-i n 

i^o) o_n 


y( 28 + n + 2,-1M 3 k 3 aS 3 ) 
2l+n„ 


(M 3 k 3> JUO < 2M 3 ) £! (n + «•) 

Upstream Shock Quantities 


The unsteady pressures and velocities ahead of the shock wave can be 
obtained from reference 1 as 
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v”* ■ E 

n=0 


1 (M,k -X )x 1 (H,k.+x )x _ 

i 1 1 n R* + e R„ 

n n 


cos 




u^e 


.lot 


E 

n=0 L 


'< M l k l - V x „+ . '<WV X - 


q n * e 


-<4 [ H o tA o< x - d o>] f l <y;H l“l> 

- 1A o[i^“? f l< yiH l“l> * “l f 2< yjH l“l>] 

cos ( “*> ♦ 1A 0 “)f 1 (»; M i“i) 


'i ■ 57 [pi (x - y) * u i (x - y) J 


In order to bring these quantities Into a form appropria te for th e shock con- 

V 2 3 

a - l^2 ( s 3 “ ¥3) 

Carrier's solution and transform them to Bessel series. 


On the shock locus 


x = (s - y)cot x 

y = ^ 


an argument of the functions f-| and f 2 of $b becomes 

i"l“l< y - s > ’ q f k 3 - 63 (s 3 - y 3 ) * <<IC3M3<S3 - y 3 > 

where 


M u 

± -7— 1 sin x + Mo cos(x - f>) 
+ k 3 j 


q l = 


..2 _ i _2 


M3 sin (x - «) 


If the generating function for the Bessel functions, and Its parametric 
derivative, 


-IMn 


w 

= J 0 (n) + ^ |t v + (-t)" v |J v (n), t = -1 (M + B), B = Vm 2 - 1 

v=l 


Ine 


+1Mn 


00 

V =1 


t = 1(M + B) 
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are Introduced, It follows that 


i ,M l“l<y- s > '“ k 3 H 3 <s 3- y 3> 
e = e 


CO 

v« ♦£ [ (t ± >B * 

m=l L J 


with 


t ± = ^ iV, 2 -1 j 


and that 


Ik. 


Vl - sin 2 ( x - 5) + 1 M,u, ( y-s) 

~ xe ' 


cos x 


1 



1 “ k 3 M 3 (s 3 _y 3> V 1 m |> 

6 L fi 


m=l 


- (-t ± )- m ] J B (« 


Consequently, the functions f-| , f 2 , appearing In $ 5 , can be transformed to 
Bessel series by using the Identities 


cos H^u^(y - s) = e 


1 ak 3 M 3 (s 3 y 3 ) 


E 

m=0 


(1 Vj m 

m 


sin M 1 u 1 (y - s) = 


1ak 3 M 3 (s 3 -y 3 ) lV^(l) 

6 lL 

m=0 


e 


m 




cos M-j co-j y = e 


1ak 3 M 3 ( s 3 y 3 ) 


E 

m=0 


( 1 *e + J ( 5 ) 

m m' 9 ' 


sin M-ji^y = - 1 e 


iak 3 H 3^ S 3 _y 3^ 


E 

m=0 


(1) ej m (5) 

mm 


x cos M-j 03-j ( y - s) 


= e 


1ak 3 M 3^ S 3 -y 3^ 


E <1> ' ! '» j »(5) 


m m' 


m=l 
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The quantities ^ are defined In appendix B.^ 

m m m m 

The Fourier cosine series appearing In each of p-] , u-| , and v-] Is trans- 
formed by again Introducing coordinates of region 3 to obtain 

[ M l k l + X n *] x " = r n 5 + ak 3 M 3 (s 3 " y 3 ) 


Then, for example. In the pressure the Fourier series becomes 



We approximate the summation In this last form by a least-squares fit to a 
polynomial of order 50, say, 


2 The functions f-j , f 2 and expressions Involving them are also listed In 
appendix B. 
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50 

2>. ,) (!' 


m=0 


m 


The polynomial Is, In turn, converted to a Bessel series by using Gegenbauer's 
expansion (ref. 10) 


(1 2 y . ^ » 2 n)r ( , * n) 


n=0 

with the result for the Fourier series In the pressure p, 


m=0 L 


1(M,k,-\ )x . 1(M,k,+\ )x 1 / \ 

• 11 " 11 " 


[n/2] 


E4 

n=0 m=0 


(1 ) nr(n - m) 
2m m! 


The alternate expression for the upstream pressure for application In the shock 
relations Is, then, 

1 [n/2] 

E (l) nr(n - m) 

^n-2 
m=0 


lut 1u 3 al V S 3 y 3^ /R 3 


P-,e = e 


E v« 


n=0 


-2m m! 





■:] 


_ 1ak 3 M 3^ S 3" y 3^ 

- e 


Ev« <1) 4 1) 


( 10 ) 


n=0 


and similarly for the velocity fluctuations 
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u^e 


io>t 


= e 


1w 3 aM 3^ s 3~ y 3^®3 


E«. 


u) 


n=0 


+ 1 A 


[n/2] 

E 

m=0 


(2) nr(n - m) 
^n-2m m! 


e W - ‘V 

n n 

o H-j slnCM^w^s) 


_ 1ak 3 H 3^ s 3 _y 3^ 
= e 


00 

E^ 

n=0 


<W (1) ^ 2) 


(11) 


/[n/2] 


V 


lut 




n=0 \ m=0 


,( 3) nr(n - m) 
2m m! 


u, ( H -Ad) 

° 0 0 ♦ i*„ i 1 

0 du 


slnfM^u^s) 


_ [ "I — 1 [ 
1 JjIntM-jCo-, s)J L 


ie w (1 >e- - (’> e . 




A 

o 1 


s 1 n( M-j co^ s ) 


-M, ^ [.1- <’>?;. <’Vj 


"V-s 0, «; 

5 e ^3"3 (s 3-»3> 0)^(3) 


( 12 ) 


Recursion Relations 

To illustrate the method of preparing the downstream quantities for sub- 
stitution into the shock relations (2), (3), (4) we compute the pressure along 
the oblique shock wave in its proper form. 
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From equations (5) and (9) 


3 ' s 


?- - 


1 k 3 H 3^ X 3 + 0 tS 3^ 


3 ax. 


co 

n=l 


cosh ne J n (k 3 r) ♦ - 1 a > 3 - ^ # t 


By employing the relations 


ar 

3X 3 


cosh e, 


ae _ slnh e 
3x 3 * - r 


Vi <x) * 'p.1 1 " 1 -f V x) 

and the derivative of the Bessel function we obtain for the derivative In the 
first term of P 3 , 



1k 3 M 3 (x 3 ,aS 3 ) 



cosh ne J (k„r) 
n 3 


IkgMgCXs+aSg) 


X 



cosh ne J n (k 3 r) 



J n 1 ( k 3 r ) cos h(n - 1 ) e 


- 2^2*" J n*l (k 3 r > cosh(n * 1)6 
n=l 


Introducing a new summation Index N = n - 1 In the second term In the brace, 
N = n + 1 In the third term, and adding the time-derivative term yields 


- Or* 


3 \ 1k 3 M 3^ X 3 +aS 3^ 


^ ^A n cosh ne J„(k,r) = ^ e 


kg _ 1 k 3 H 3 (X 3+ aSg) 
n'^' ' " 2 


n=l 


OO 

^ ] cosh ne J n (k 3 r) 


n=0 


Vl (1 ‘ 4 n0 )(1 “ 6 nl } ~ Mg A n (1 " W “ A n+1 


For the term Involving §b we again use the generating function of the 
Bessel functions to obtain 
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- 1 f W y 3 -S 3^ _ 1 ak 3 M 3^ S 3 y 3^ 1 q 3 5 


1ak 3 M 3^ S 3“ y 3^ 


J o<« + 


GO p 

( t-) m * ( -ti; 
m=l ■ 


tV m 3.(0 


and, on the shock wave, 
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From these expressions all of the functions of X 3 and y 3 In the potential 
$b may be represented as expansions In Bessel functions In the same way as the 
corresponding quantities In the upstream region. After considerable manipula- 
tion one can write 


s_/. _3_\ - s 0 1 aM 3 k 3 (s 3" y 3 ) (3) (1). , c . 

=3( 3 - ax 3 ) * b -ip Zf ** V« 


with ^ ^defined as In appendix C. Thus, on the oblique shock where 


x 3 + ttS 3 = a(s 3 - y 3>’ 
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Likewise, the unsteady velocities and shock displacement are given by 
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We can now substitute the expressions for p-j , u-| , P 3 , U 3 , V 3 Into the 
shock relations (2), (3), (4), multiply the first equation by 2 s 3 M 3 /sk 3 , the 
other two by 2 /k 3 , and by equating coefficients of the Bessel functions of 

like order obtain the following recursion relations for the A . C . and D : 

n n n 


where 
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are defined In appendix C. 


The matrices 
The recursion 


^y and the vector 3 ^ 

system Is entered at n = 0 with 


= C -1 - D -1 


- 0 - A 0 - C 0 


and with D q given by equation (8). 


Channel Region 

The wedge-flow solution Is valid up to the point where the shock wave In- 
tersects the adjacent blade. At the Intersection the flow Is at the angle 63 
(see fig. 3 ). We construct a perpendicular to the blade at this point (In the 
y3-d1rect1on) and designate the region downstream of It up to the normal shock 
wave as the "channel region". The passage diverges slowly from Its Initial 
width S(X3) = S3 to Its final value S(X3) = S2- Provided dS(x3)/dX3 « 1 , 
It Is appropriate to represent the mean flow as one-dimensional, l.e., 

^3 =^3(X3). For the unsteady flow also to be slowly varying It Is nec- 
essary that the wavelength of the disturbance X ~ S(X3). 

The slowly varying approximation for the unsteady flow Is a straight 
channel solution modified to account for the slowly varying channel width. The 
potential and stream function respectively, are 


*3 " 


CO 

E 


n=0 


r 


T n e 


3 ( M 3 k 3-X< 3) )d x 3 l/' X3 (M 3 k3 + x( 3 ))dx 3 

+ T n e 


cos [—) 


+ *b + *s 


( 14 ) 


¥ 


3 





( 15 ) 


where the parameters M 3> g> 3 = uc/# 3 , k 3 = “3^/83. and 


4 3) * 


( mr/B^S) ^ vary slowly with x 3> 


The additional potential $ s , appearing In equation ( 14 ), arises from 
the motion of the upper blade surface In following the disturbance of the 
oblique shock wave. In figure 5 the dashed line represents the displaced 
position of the model wall due to a displacement x s (0) of the oblique shock 
wave from Its mean position. The displacement of the wall In region 3 normal 
to Itself Is 
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where from equation (7), 


W 


s 


x $ ( 0 ) sin 5 
sin x 


lak ( 0 )M ( 0 )s„ * 
x $ (0) = e J d J D 



The upwash condition on the blade, 


( 16 ) 


Is satisfied, along with the convective wave equation governing the flow 
(eq. ( 2 . 2 ), ref. 1 ) by 

cos[M 3 ( 0 )co 3 ( 0 )(s 3 - y 3 )] 

$ s = " 1 u l W s o> 3 ( 0 )H 3 ( 0 ) s 1 n[M 3 ( 0 )« 3 ( 0 )s 3 ] 


+ (2) 

In order to find the expansion constants T" , ' In equations (14), 

(15) we match the latter quantities to the wedge solution (5), ( 6 ) at X 3 = 0. 
To expedite this we can, without affecting P3, U3, V3, replace $3 and ¥3 
In the wedge region by 
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~3 

The vortlclty, , contributes no normal velocity at the blade surface, a 
characteristic shared by equation (15). 
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In order to match pressures at X 3 = 0, It Is necessary to match both §3 
and 35/3x3 . Equating §3 and $3 at X3 = 0, multiplying by sin miry3/S3, 
Integrating and employing orthogonality properties of the trigonometric func- 
tions yields 


T_+ T~= 
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The Integral A n>m appearing In In Is evaluated In appendix 0. 
Similarly, matching 35 / 3 x 3 and 3 $/ 9 X 3 yields 
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In the second term of the Integrand, and 
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In the third term, and 
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Introducing a new summation Index N-j = n - 1 In the second term of the Inte- 
grand and N 2 = n + 1 In the third term leads to 
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When equations (16) and (17) are combined one finds 


e 1H 3 (0)k 3 (0)s 3 OL ® 

n = s 3 (1 + Vo* Z 


1k 3 (0) 


^ A m 1 2 \( 3 ) ( 0 ) (Am+1 " * m ’ n 


2Xj u '(0) 


+ 11 + 


k 3 (°) 


o 3 (0)C 


M 3 (0 )X n 3> (0) / “3 (0) + < nir/s 3 > 2 


7 ¥ 0)k 3 (0) \ M 1 s1n 6 D 

\ 1 x n 3) <°> / ( nir/s 3 ) 2 - M 3 (0)^(0) sin x 


22 



Equating vortlcltles ^ and at x 3 = 0, multiplying by 

sin miry 3 /s 3 = (-l) 111 ^ s1n(s 3 - y 3 )mir/s 3 and Integrating across the channel 

yields 
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Is evaluated In appendix D. 

The pressure obtained from the potential $3 of equation (14) Is, on the 
upper surface of the blade. 
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with x 3 = x 3 - d s + s 3 - In calculating the lifts from these pressures we must 

allow for the difference In the slope of the curved blade surface In region 3 
and the slope of the blade surface In region 1. On the upper surface the 
contribution to the lift Is 
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Defining the moment with respect to the trailing edge, we have 
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On the lower surface of the blade we also allow for the altered length of 

the channel (d a versus d s ) and obtain for the lift 
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The corresponding moment about the end of the channel Is 
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Region 2 

The solution for the unsteady flow In the subsonic region that follows 
the normal shock wave proceeds as In reference 1 after allowing for the down- 
stream displacement of the shock to the trailing edge of the blade and Includ- 
ing the effects of upstream vortlclty. The downstream pressure Is still given 
by equations (61) and (G2) of reference 1, but now the coefficients Bf, 2 ) of a 
downstream propagating wave are determined by the equation 


s 2 1 


+ 6. 


B< 2 >e 




+ a. 


Vo 


n,o 


GO 
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m=0 


B* 2) K 
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mn 


rather than by equation (3.33) of reference 1. The small quantity e 0 repre- 
sents the length of a nominal channel between the shock position and the end of 

+ 

the blade. The quantities a. K remain as defined In reference 1, and F 

n mn n 

Is defined In appendix E. 

The pressure on the blade surface In region 2 Is given by equations (Gl) 
and (G2) of reference 1. The lift and moment derived from them are Included 
In appendix E. 


Shock Lifts and Moments 

The unsteady displacements of the two shock waves where they Impinge upon 
the blade surfaces contribute to the unsteady lifts and moments by shifting the 
boundary between adjacent regions of the flow and thereby altering the domain 
of Influence of the mean pressures, as discussed In reference 1. The oblique 
shock, for example, forms the boundary between regions 1 and 3 on the top of 
the blade. From figure 5 It can be seen that x s (0)/s1n x Is the penetration 

of region 1 Into region 3, for which (pj^ - p 3 °^ )x $ (0)/s1n x Is a correc- 
tion to the downward pressure force normal to the surface In region 1. Since 
the dimensionless pressure difference Is 


26 




It follows from equation (18) that the lift due to the oblique shock Is 



sin <5 Ds 
sin x cxS^ 


1aM 3 (0)k 3 (0)S 3 


The corresponding moment Is 

= (s cot x - V-^s^ 

where d 0 Is the value of x at which the elastic axis projects upon the sur- 
face In region 1 . 


The expression for the lift due to the normal shock Is similar to that 
given In reference 1 with modifications required for the downstream displace- 
ment of the shock and Its Impingement on the lower surface rather than the 
upper: 
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The moment Is 


^ 2) - [m' - f > c ° s < b 2 - *»> - .; K : 


RESULTS AND DISCUSSION 


A set of six blades which are examples of the two-shock model Is shown In 
figure 6. For all blades the solidity Is 1.4, the stagger angle Is 62° and 
the entrance angle Is 67°. The exit angle Is 65° except In the first case for 
which It Is 63°. The Initial choices of entrance Mach number, M-| , and pressure 
ratios, P 2 i , were made to provide a distribution of flow conditions over the 
shaded region of figure 4. Since figure 4 was constructed for the simplest 
geometry having straight surfaces and a boat-tall trailing edge (dashed lines 
In fig. 6) It Is necessary to recompute the pressure ratio after the choice of 
boundary-layer thickness and exit angle to account for the Isentroplc expansion 
of the gas In the channel and the strengthened normal shock wave. The expan- 
sion and the enhanced normal shock strength tend to counteract each other; thus 
the actual pressure ratios stay within 15 percent of the Initial estimate and 
are sometimes much closer. 
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Each blade has been given a traillng-edge thickness to represent the dis- 
placement effect of surface boundary layers. The displacement thickness on 
the upper surface Is denoted by <$ u 2|j and the displacement thickness on the 
lower surface by The values for S u and were chosen to yield a 

combined thickness of about 10 percent of lb- In one case 6 U = and In 
two cases 5^ = 0. In the remaining three blades Sj Is actually negative, 
thereby reducing blade thickness, but simultaneously Increasing the effective 
stagger angle. This suggests that In some applications It may be desirable to 
begin the model design with an Initially undervalued stagger angle and to use 
negative boundary-layer thickness and the exit angle to bring It up to the 
desired value. 

Some representative pressure distributions on the blade surfaces are shown 
In figures 7 to 12, for blades pitching about their centers (f = 0.5). The 
pressures for four Interblade phase angles are given In each case. On the 
upper surface, the pressures In the entrance region exhibit the same charac- 
teristics that occur In the normal-shock configuration of reference 1, namely, 
that they are relatively uniform and low In magnitude over the full length of 
the surface -s f < x < s cot x except, perhaps, for phase angles of ir/2 or ir, 
or at high frequency. The pressures on the additional surface 0 < x < s cot x 
now Included In the entrance region have no unusual characteristics. In the 
second supersonic region behind the oblique shock wave, however, the pressures 
are sharply peaked, a characteristics observed previously (ref. 1) In the sub- 
sonic region behind a normal shock wave, and now seen to be feature of the un- 
steady flow behind a weak shock as well. These peaks appear to be even more 
predominant at low frequencies than at high. 

On the lower surface, as well, the unsteady pressure behind the oblique 
shock wave has generally several times the magnitude of the pressure In the 
entrance region, although It Is likely to be more uniform In spatial distribu- 
tion except at the high frequency. At the point on the blade surface where 
the curved surface of the channel begins, the pressures are usually strongly 
affected. There Is often a sharp reversal In the magnitude of the pressure 
which then falls to a lower magnitude at the end of the channel, just preceding 
the normal shock wave. 

When the disturbances are transmitted through the second, normal, shock 
wave there Is once again apt to be an amplification of their strength, espe- 
cially when the blades are out of phase (a = ir/2) or of opposing phase (a = ir) . 
When the blades are moving In phase (a = 0) the pressure amplitude In the sub- 
sonic region Is likely to be low. In general, however, the effect of succes- 
sive shock waves appears to compound, thus raising the disturbance magnitude 
to a high level . 

The dimensionless lifts and moments In regions 2 and 3 have been derived 
or recorded In the Analysis section. The lift and moment In region 1 can be 
obtained from the pressures given In equations (Dl) and (D2) of reference 1. 

They are defined and recorded here In appendix F. The separate lift forces 
when added vectorlally and referred to the normal of the entrance-region sur- 
face yield for the total lift 

1 1 m ( 2 ) 1 1 

L = (l_ 33 + L 32 - L 32 ) cos S + l 21 cos(B 1 - I3 2 ) - L n + L s + L s cos ( fi i ~ ^ 
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Likewise, after each moment M^ has been translated to a corresponding moment 
referred to the elastic axis, the complete moment Is given by 


M = M 00 + M„„ + M 


'33 


32 


21 






The total complex moment M Is plotted for each of the blade models at three 
reduced frequencies, «■] , In figures 13 to 18. When 9m M In these figures 
Is negative, the blade is receiving energy from the flow and flutter Is the 
result. The Interblade phase angle a Is the parameter along the curves. 

Black dots denote an Interval of ir/8 In a. 

In all cases the figures Indicate that there Is Instability In torsion at 
the moderate and low frequencies but In only one case at the high frequency. 

The exception Is blade 4 which has the largest thickness ratio of all the 
blades. This progression to greater Instability as the frequency Is decreased 
was also the characteristic of the cascade with one normal shock wave In 
reference 1. The Instability Is more apparent In the double shock case with 
flow at the Intermediate frequency being decisively unstable. The Increased 
complication of the flow In the double-shock model Is reflected In the more 
complicated stability figure than for the single-shock flow. The number of 
loops and crossings has Increased and the moment wanders more erratically In 
the complex plane. 

The complex lifts are displayed In figures 19 to 24. A negative 9m L 
means that the blade Is receiving energy from the flow and there Is Instability 
In bending. Qualitatively, the trend towards Instability In bending at lower 
frequencies Is the same as In torsion. Blade 4 Is seen to be the most stable 
of the six examples, and, again, the figures display considerably more struc- 
ture than In the single-shock model. 


CONCLUDING REMARKS 

The model for a supersonic blade row with shock waves which has been pre- 
sented here has extended the operating range of turbomachinery for which 
stability calculations can be made. Based on the single geometry explored 
here, the model Is appropriate at M-| >1.3 at the lowest applicable stage 
pressure ratio and, for any pressure ratio, at higher entrance Mach numbers 
than the single shock model. The calculations for specific cases show that In 
this operating region the presence of shock waves has a strong Influence on 
the unsteady flow and that the tendency for Instability at low frequencies 
observed In the single-shock model of reference 1 persists In the two-shock 
model. Thickness ratio appears to have significant Influence on stability of 
the flow at high frequencies. 
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APPENDIX A - COEFFICIENTS OF SHOCK CONDITIONS 
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APPENDIX B - FUNCTIONS IN THE UPSTREAM SHOCK QUANTITIES 
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APPENDIX C - MATRICES OF THE RECURSION RELATIONS 
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APPENDIX D - INTEGRALS ARISING IN THE MATCHING PROCESS 


In the integration across the blade passage at X 3 = 0 there arise 
Integrals of the type 
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puts these Integrals in the form 



In succession one expands the Bessel function in its power series, translates 
to a new coordinate rj = n ± X and converts factors ( rj + 2\) K to their 
binomial expansion (\ = b/p). Then one has Integrals of a standard type 
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To put A mf0 In a form convenient for calculations we first Introduce a 
new summation Index 2. = k - j and sum over all nonnegative values of both a 
and j. In the a,j matrix we construct a final Index k = 2a + j which sums 
over all nonnegative values with 0 < a < [ k/2 ] . Then 


m 


a 3 y 

,n " 2nir 1 — 4 
j=0L 


0>*> * 0--> 


A^ m) (b) 


r< n) (v) 




ki (-1)' 


2r+l 


(k - 2r - 1) 


(1 + v) 


k-2r-l 


Aj m) (b) 


[J/2] 

- £ 

k=0 


(b/2) m+2 ^~ 2k 2 j_2k 
(m + j - k) ! k! (j - 2k)! 


For the case n = 0, 


r k (±y ’ ntr) _ (1 ± u ) k+1 
nw ± v( k + 1 ) 

Another quantity arising from the Integration across the blade passage Is 
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y Is the Incomplete gamma function. 
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APPENDIX E - EFFECTS OF UPSTREAM VORTICITY AND CHANGED 
SHOCK LOCATION ON SOLUTION IN SUBSONIC REGION 


Equation (3.33) of reference 1 must be modified to account for the 
presence of vortlclty upstream of the normal shock wave and for the altered 
position of the shock wave. 

If vortlclty ijr] Is present. It affects the derivation of the shock con- 
dition (A2) (or (2.18)) as given In appendix A of reference 1. The expressions for 
v-| and 3<}>-|/3t must be augmented to read 
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From these one can form a(v 2 - v^/at, employ D^/Dt = 0 on both sides of the 

shock wave, and set the expression equal to the time derivative of equation 
(Ale) to obtain. 


a_ 

ay 


*2 /' 
p 2 +U 2- 


aX i 

3t 


1 3X S 

+ p i + u i " iF 


°u 2 /3> 2 a> 2 \ r *1 3 V 


<%- 


i \ay 


3t 2 


3t 


When (Ala) and (Alb) are Introduced, the augmented version of equation (A2) Is 
found : 
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When the time differentiations are carried out and the time factor Is can- 
celled, we have an equation corresponding to equation (2.18) of reference 1, 
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Again In the companion shock relations (eq. (2.17) of ref. 1) we allow for the 
occurrence of upstream vortlclty and write 
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If ¥ 2 Is eliminated between (E2) and (E3) by differentiating (E2) with 
respect to y and obtaining the first and third y derivatives of Yp 
from ( E3) , the extra term In In the combined equation Is 
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we find the vortlclty term to be 
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In reference 1 the region before the normal shock wave was region 1; here 
It Is region 3, and we substitute for V-\ the expression In equation (15) for 
¥ 3 , evaluated at X 3 = d s : 
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Substituting M_(d ) for M, everywhere In the shock condition and Introduc- 
Ing the shock position x = 1 - s 2 - c Q Into the expressions for the poten- 
tials where they appear In the shock relation, we find the condition on the 
coefficients B^ of the duct waves In the subsonic region. 
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Equation (E4) replaces equation (3.33) of reference 1. In It e 0 Is the 
length of a nominal duct ahead of the trailing edge In which the flow Is sub- 
sonic. In the calculations made for this report e 0 has been set to 0.01. 
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There Is a factor s/S2 appearing In L21 and M2i which arises from 
the normalization of the pressure In region 2 and which merits explanation. 
Because time, velocities and pressure are nondlmenslonalized with respect to 
the chord, the density p, and velocity °U-\, whereas the material derivative In 
region 2 Introduces the velocity #2* the ratio of mass fluxes p^/pi^ 
appears In the nondlmenslonalizatlon of the momentum equation. Where the blade 
passage Is of constant width, as Is the case In reference 1 , the ratio Is 
unity; but when the blade passage varies In width, as It does here, the ratio 
Is just the reciprocal of the channel widths, s/S2- This explains the appear- 
ance of the latter ratio In the lift and moment as well as the appearance of 
S/S3 In equation ( 9 ). 
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APPENDIX F - LIFT AND MOMENT IN REGION 1 
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Schock location determined 

from the LA measurements, 
15 percent span 

Schock location determined 

from holography, 5 percent 



Figure 1. - Relative Mach number contours measured with the LA at 
15 percent span from the tip (ref, 3). 



Figure 2. - Rotor blade tip static pressure contours. Inlet Mach number. 1. 493 
(ref. 5). 
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Figure 5. - Displacement of blade wall due to incremental change 
in oblique shock position. 






estimated P 21 - 2- 75- 



Figure 6. - Blade models for six flow conditions. Entrance angle, 67°; stagger angle, 62°; solidity, 1.4. 
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d s -l -S -S3 +d i 0 s 2' s 3 +d i 

(d) Out-of-phase component (imaginary part) of pressure acting on lower surface. 

Figure 7. - Concluded. 
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Figure 8. - Nondimensional surface-pressure amplitudes on blade 2 undergoing pure 
pitching about its center; reduced frequency (l/2)«j ■ 0.5. 
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(c) In-phase component (real part) of pressure acting on lower surface. 
Figure 8 . -Continued. 
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Figure 9.- Continued. 
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(a) In-phase component (real part! of pressure acting on upper surface. 

Figure 10. - Nondimensional surface-pressure amplitudes on blade 3 undergoing pure 
pitching about its center; reduced frequency (1/2) wi - 1.0. 





(b) Out-of-phase component (imaginary part) of pressure acting on lower surface. 
Figure 10. - Continued. 
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(d) Out-of-phase component (imaginary pai 
Figure 10. - C 
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Figure 11.- Continued. 
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(c) In-phase component (real part) of pressure acting on lower surface. 


Figure 11. -Continued. 
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(a) In-phase components (real part) of pressure acting on upper surface. 

Figure 12. - Nondimensional surface-pressure amplitudes on blade 5 undergoing pure 
pitching about its center; reduced frequency (1/2) Uj - 0.5. 
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(d) Out-of -phase component (imaginary part) of pressure acting on lower surface. 
Figure 12. - Concluded. 
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-2 0 2 4 

10 Uj - 2. 

Figure 14. - Moment coefficient for pitching motion about center of blade 2. (Mj ■ 1.6; P 2 i' 2.76.) 
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(c)idj- Z 


Figure 15. - Moment coefficient for pitching motion about center of blade 3. (Mj • 1.6; P 21 ■ 2.81. ) 
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(clUj- 2. 

n about center of blade 4 (Mj-1.8; p. 


- 3 . 49 .) 
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(c) wj ■ 2. 

Figure 20. - Lift coefficient for plunging motion of blade 2. (Mj • 1.6; 
P 21 - 2.76.) 
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